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Tin's work is the first step toward increasing the bandwidth of the DSS-13 an- 
tenna position loop controller. A wider bandwidth of the controller allows for faster 
tracking rates and better pointing performance under windy conditions. To achieve 
this goal, the antenna control system model has to be improved, such that it will 
accurately reflect the dynamic properties of the antenna. The existing analytical 
model, due to its many uncertainties, could not be used in the design of the con- 
troller. However, by using experimental data, the analytical model is modified and 
improved, and a new model is obtained through system identification techniques. 


I. Introduction 

The position controller at the DSS-13 beam- waveguide 
antenna does not employ model-based control law. In the 
interest of implementing model-based control law at the 
DSS-13 antenna, an analytical model of this antenna was 
developed [1]. The experiment described in this article 
was then designed to validate this analytical model and 
to identify a new antenna model from the experimental 
data (the latter model will be referred to as the experi- 
mental model). The data collected by this experiment at 
the DSS-13 antenna were collected for each axis separately. 
The data were processed and used to modify the analytical 
model and to obtain an experimental model through sys- 
tem identification techniques. Since the latter approach 
requires only a portion of the data, a sample was taken 
from the remaining data and used to validate the experi- 
mental model. 

II. Description of the Experimental Software 
and Hardware 

The purpose of this experiment is to gather time series 
data on the input to the antenna drives and the output 


of the antenna position. These data are processed to de- 
termine the frequency responses, or transfer functions, of 
the DSS-13 antenna. These transfer functions represent 
an open-loop system. They are used as a reference to de- 
termine the correction of the parameters of the analytical 
model, and also to identify the transfer function of the 
system using identification software. 

A. Experimental Software 

The LabView 2 application software generates an ana- 
log random (white noise) signal. This signal is input into 
the rate loop via the test input lead in the servo interface 
chassis, and the encoder output is sampled from the data 
converter. The input and the output signals are buffered 
until the designated number of samples is taken. The sam- 
pled data are then saved in floating point MATLAB 1 load- 
able format files. 

B. Experimental Hardware 

The data acquisition system consists of the LabView 2 
software running on a Macintosh Ilfx computer. The com- 

1 Copyright 1985-91, Inc. All rights reserved. 
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puter contains an analog board which can output analog 
signals through two channels to the antenna rate loop. It 
also contains two digital boards which sample the azimuth 
and the elevation position encoder output (Fig. 1). 

The National Instrument NB-MIO-16 multifunction 
analog board for the Macintosh is used to produce the 
analog rate signals injected into the antenna rate loop. 
The board contains two 12-bit digital-to-analog convert- 
ers with voltage outputs. These two analog outputs are 
hooked up to the test input leads of either the azimuth 
or the elevation rate input, or both. The test inputs are 
located in the servo’s interface chassis. The servo inter- 
face is switched to manual mode so that a rate signal is 
accepted from the test input lead instead of the antenna 
position controller. The brakes in both axes are released 
manually before the experiment begins. 

The National Instrument’s NB-DIO-32F 32-bit, paral- 
lel, digital, input/output interface board for the Macin- 
tosh II is used to sample the binary output of the data 
converter. The data converter buffers the 24 bits of en- 
coder data available in each axis and provides the data in 
latched binary form. 

A Hewlett Packard frequency counter is used to mea- 
sure the sample period for each run of the experiment, 
since the sample period depends on the execution speed of 
the LabView 2 software, and it cannot be determined in 
advance. 

C. Data Acquisition 

The experimental data were collected separately for 
each axis, namely, for elevation rate input to elevation 
position output, elevation rate input to azimuth position 
output, azimuth rate input to azimuth position output, 
and azimuth rate input to elevation position output. The 
maximum rate analog input voltage was 1.3 V, and the 
minimum was 0.1 V. 


and 3.6 m/sec). When no input was applied to the eleva- 
tion axis, data were taken for the elevation fixed at both 
90 and 60 deg. When no input was applied to the azimuth 
axis, it was left fixed at 140 deg. 


III. Analysis 

A model-based control is currently not employed at 
the DSS-13 antenna. In order to design a model-based 
controller, a suitably accurate antenna model is required. 
There were two approaches taken to obtain such a model. 
First, the model presented in [1] was adjusted to fit the 
data collected in this experiment. In a second approach, a 
model was obtained by using a system identification tech- 
nique based on the data. 

A. Model Description 

There are two analytical models for the DSS-13 antenna 
[1,2]: a model for the antenna at 90 deg elevation, and the 
other for the antenna at 60 deg elevation. These models 
are described in the state space form 


x = + y = Cx (1) 


In this model, the state vector x is of dimension n, the 
input u is of dimension p, the output y is of dimension q ) 
and the matrices A , £, and C are of dimensions n x n, 
n x p, and gxn, respectively. The full order of the analyt- 
ical model in [1] is n = 90, and the reduced order model 
is n = 27. For the system whose input is azimuth and 
elevation rate and whose output is azimuth and elevation 
position, p = 2 and q = 2, and 


B — [B az B e i], 



( 2 ) 


For each axis and each of its cross couplings, the ex- 
perimental data were taken in separate sequences. Each 
sequence contained 8192 data points, and there were any- 
where between 10 and 40 sequences per run of the exper- 
iment per axis. The available travel in the elevation axis, 
in conjunction with the rate input and the sample interval, 
restricted the number of sequences that could be taken per 
run of the experiment. The azimuth axis, with 450 deg of 
travel from one limit to the other, never proved to be a 
limiting factor given the same rate input and sample pe- 
riod. 


where B az and B e i represent azimuth and elevation input 
rate, and C az and C c i represent azimuth and elevation 
encoder position. 

The four transfer functions (elevation rate to elevation 
encoder position, elevation rate to azimuth encoder posi- 
tion, azimuth rate to azimuth encoder position, and az- 
imuth rate to elevation encoder position) are used to com- 
pare analytical and experimental data. Let Tij(s) be the 
transfer function defined by 


The data were gathered on several occasions. On each 
occasion, the wind speeds at DSS 13 were low (between 0 


Tij(s) = Ci(sl- A) l Bj, i,j = az,el (3) 
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By judiciously varying a set of A, and Cj matrices, one 
changes the transfer functions such that they better fit the 
empirical transfer functions. A “good” fit is required only 
over the frequency range where high coherence between in- 
put and output is observed, in most cases for frequencies 
from 0.01 to 10 Hz. In this case, the good fit means that 
the peak frequencies and magnitudes of the model trans- 
fer function line up with those of the empirical transfer 
function. 

B. Frequency Responses of the DSS-13 Antenna 

Using frequency domain analysis, the data were used 
to modify the analytical model. The input, u(t), and out- 
put, y[t ), time series data have been detrended and passed 
through a Hanning filter to prevent spectral leakage [3]. 
Then, a fast Fourier transform (FFT) was performed on 
the resulting time series data. The magnitude of the trans- 
fer function, T' uy (/), and the coherence, / (/), were esti- 
mated from the time series data (c.f., [3]) using u(t) and 
y(t), filtered and detrended input and output vectors of 
8192 samples each. 


Tuy(f) = 


Puyjf) 

Puu(f) 


(4) 


^ _ 1P UV (/)1 2 
^ " Puu(f)PyyU) 


(5) 


C. Modification of the Analytical Model 

The analytical model as described in [1] is a combina- 
tion of the antenna structural model and models of the 
elevation and azimuth drives. It is derived from the best 
available knowledge of the antenna structure and servos. 
The antenna structural model is obtained from its finite 
element model. Although complex, as the finite element 
model usually is, its accuracy is still limited. For exam- 
ple, damping ratios or “non-st ructural” masses, such as 
counterweight mass, are usually only roughly estimated. 
The available field data allow one, to some extent, to cor- 
rect some parameters of the analytical model so that its 
properties fit more closely the properties derived from the 
experiments. In addition, the drive models have some un- 
certainties. The gearbox stiffness is not known precisely, 
since it depends on the countertorque value, which is itself 
a fuzzy number. Also, the gains of the drive amplifiers 
are not set precisely, or can change. These and other less 
known factors impact the model accuracy. 

D. System Identification Model 

A second model is identified directly from the experi- 
mental data using the system identification software Sys- 
tem/Observer/Controller Identification Toolbox (SOCIT), 
written in the MATLAB language [4]. The “okid” func- 
tion of SOCIT identifies a state space model, i.e., the A, 
B y and C matrices, given the input-output data, sample 
period, and the number of observer Markov parameters. 
For more detailed information about SOCIT software and 
this function, the reader is referred to [4] and [5]. 


where P U u(f) is the power spectral density estimate of 
u(t), Pyy{f) is the power spectral density estimate of y(t)y 
and P U y(f) is the cross spectral density estimate of u(t) 
and j/(<). The phase is estimated from the time series data 
by using the following formula: 

^ (/) = arctan ^(ij (6) 

where i?e(.) and 7m(.) denote real and imaginary parts. 

The transfer function estimates were obtained by aver- 
aging the magnitude and phase of each sequence. The 
magnitude, T uy {f) } and the phase, are plotted 

versus frequency in Figs. 2(a)-(f) (dashed line), where 
/ = w/2tt is frequency in Hz, and u is frequency in rad/sec. 
In Figs. 3(c), 4(c), 5(c), and 6(c), the coherence was also 
plotted to determine the range of frequencies over which 
the data were valid. 


The SOCIT software was used to identify a model for 
each axis and cross-axis, so that the A , B , and C state 
space matrices were obtained for each subsystem. The or- 
der of the system was chosen based on the system’s Ilankel 
singular values plotted as part of the okid function output. 
Typically, the order chosen was between 25 and 30. The 
resulting state space matrices were transformed into bal- 
anced coordinates, so that a matrix A was in diagonally 
dominant form. The 2x2 diagonal blocks represent the 
system modes in decreasing order of importance. The di- 
agonal elements of a 2 x 2 block represent the system 
damping, and the off-diagonal elements represent natural 
frequencies at those modes. 

The order of each subsystem was reduced to the small- 
est acceptable order, i.e., such that all the modes visible 
in the data up to 10 Hz were preserved. The order, n , 
of the T azaZ y T aze i , T e uh and T eiaz subsystems was n = 
14, 13, 10, and 11, respectively. The identified model had 
some discrepancies (with respect to the experimental data) 
which apparently could be removed. Namely, the damping 
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ratios were adjusted for the system modes which appear 
either underdamped or, more commonly, overdamped. 

These four identified subsystems are then used to sim- 
ulate time series responses, given the actual input data, 
u az or u c [, and small bias signals, u azb or u e i b) 


Vazaz — ^azaz u az (7&) 

Vazel ~ r E' a ztl u az H" '^ > elel u elb (7b) 

Vtui = T e ielU e i (7c) 

Velaz — T'azaz^azb + T ela2 u el (7d) 


The use of small bias in the cross-coupling responses was 
necessary because the straight gains | T azaz [ and | T t \ t \ | 
are much larger than the cross gains | T aze \ | and | T t \ az |. 
For the same reason, the bias can be neglected in the 
straight responses (azimuth- to- azimuth and elevation-to- 
elevation). The simulated responses, y 8 azaz , y 8 aze ^ VeieV 
and y* elaz , were compared to the actual SISO output data, 
Vazaz , Vazel, Veiei, and y e i az , and the approximation errors 
in azimuth and elevation, e azazy t aZ ei, teUh and t t \ aZy were 
computed. 


tazaz — 


Vaza 1 J/LaJI x 10Q% 


II Vaz 


Cazel - X 100% 

II Vazel II 


Celel = 11 yt ||‘ f y f? e ' '' X 100% 
II Velel || 


telaz — 


Velaz Velaz 
II Velaz || 


X 100% 


(8a) 

(8b) 

(8c) 

(8d) 


These formulas estimate the relative discrepancy between 
the measured and simulated signals. 


IV. Discussion of the Results 

A. Results of Analytical Model Adjustments 

The magnitude plots of the four transfer functions ob- 
tained from the analytical model (solid line) and from the 


experimental data (dashed line) are shown in Figs. 7(a)- 
(d). The discrepancies between the analytical and experi- 
mental transfer functions are immediately obvious. In par- 
ticular, there is a mismatch in the first resonance frequency 
in the azimuth-to-azimuth and elevation-to-elevation 
transfer functions, and there are higher resonance peaks 
in the analytical model for almost all resonance frequen- 
cies. The mismatch in the first fundamental frequencies 
is due to the underestimation of nonstructural masses in 
the structural analysis. The high peaks at the resonance 
frequencies are caused by an assumption of very low struc- 
tural damping (0.5 percent). Modifying the rigid-body 
structural masses (an increase of 45 percent in the ele- 
vation rigid body modal mass and an increase of 70 per- 
cent in the azimuth rigid-body modal mass) and increasing 
modal damping for the higher frequency modes to 5 per- 
cent improves the fit between the analytical model and 
experimental transfer function curves. Further improve- 
ment is obtained by adjusting the drive parameters. In 
the elevation and azimuth drives, the amplifier gain, de- 
noted k r in [1], is reduced by 30 percent (from 80 to 56), 
and in the azimuth drive, the gear box stiffness is reduced 
by 15 percent (from 2 x 10 7 to 1.7 x 10 7 ). The results of 
the modifications of the analytical model show a better fit 
to the experimental results, as seen in the solid-line plots of 
Figs. 2(a)-(d). The worst fit between the experimental and 
model transfer function curves occurs in the elevation in- 
put to azimuth output. This is due to the very small value 
of this function. It is at least 100 times smaller, in mag- 
nitude, than the elevation-to-elevation transfer function, 
or at least 10 times smaller than the azimuth-to-elevation 
transfer function. The complexity of the analytical model 
made it difficult to determine what other parameters were 
responsible for the remaining differences. 

The significant modal frequencies as estimated from the 
data are shown in Table 1. 

B. Results of System Identification 

In Figs. 3, 4, 5, and 6, the identified transfer function 
plots are presented with the empirical transfer function 
plots. Both plots show low frequency anomalies. The 
cross-coupling transfer functions T e j az and T aze j are ex- 
pected to approach zero value at zero frequency, since no 
static coupling between elevation and azimuth or between 
azimuth and elevation is observed. However, in Figs. 4(a) 
and 5(a), such a tendency is not present. The magnitude 
of the transfer function from azimuth to azimuth or from 
elevation to elevation should roll off at a 20-dB/dec rate 
for low frequencies, since the system contains two integra- 
tors (or two poles at zero). This can be observed in the 
measured frequency range of these plots. However, for fre- 
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quencies lower than measured (/ < 0.01 Hz), the identified 
transfer function does not rise. 

The departure of the empirical transfer function plots 
from the expected zero value in the low frequencies is 
caused by the presence of a small bias at the input port 
where no signal is applied. For example, when an eleva- 
tion rate signal is applied and the azimuth position is sam- 
pled, a small azimuth rate input still exists, and it shows 
up in the output data at low frequencies. Its presence is 
explained with the comparatively high gain for straight 
connection (azimuth-to-azimuth or elevation-to-elevation) 
when compared with the cross-connection (elevation-to- 
azimuth or azimuth-to-elevation). The unmeasured bias 
can be explained by the lack of causality between input 
and output in cross-coupling for low frequencies, as visible 
in the coherence plots of Figs. 4(c) and 5(c), where the 
coherence is almost zero. In order to correct this problem, 
the data were filtered with a high-pass filter before being 
used in the system identification software. As a result, the 
cross-coupling transfer function plots are obtained as seen 
in Figs. 4(b) and 5(b). 

For the straight transfer functions (azimuth-to-azimuth 
or elevation-to-elevation), the discrepancies between the 
identified model and the empirical transfer function in the 
very low frequency range are due to the presence of two 
very small eigenvalues of the identified A. By setting them 
to zero, the discrepancies were eliminated. 

The measurements were used to obtain both the ad- 
justed model and the identified model. However, the iden- 
tification software requires a relatively small portion of 
the data gathered. Hence, the remaining time series data, 
Vazaz j Vazti) Veieii &nd yeiaz> taken at the DSS-13 antenna 
are compared to the simulated time series data, y * azaz , 
vLet’ y’eleh and viiaz > of the identified model (Figs. 8-11). 
The discrepancy between the two signals, relative to the 
scale of the original signal, are computed using Eqs. (8a)- 
(8d). 

An additional comparison was performed as follows. 
Define a positive function ri(k) 

W* i = e/ - « ( Q ) 


where superscripts m and s denote measured and simu- 
lated output, respectively. For ki = 1, r< shows the dis- 
tance (or an approximation error) between the experimen- 
tal and simulated data. One can further improve the fit 
between the experimental and simulated results by varying 
the parameter ki. This additional fit is possible because 
the system was identified from a detrended set of data, but 
is compared to a non-detrended one. Let ri(ki) achieve 
the minimal value for ki = ki Q) i.e., let rj(&,*) > r(&,* 0 ). 
If k{ 0 = 1, the simulated series is the best-fitted one and 
no modification is necessary. Thus, the gain adjustment 
factor 5*, defined as 

Si =| k i0 - 1 | x 100%, i = az, el (10) 

is the measure of good fit. The gain factor 6 t * is the per- 
centage that the system gain is adjusted to improve the fit 
of the simulated time series y a zaz$ and y e j e js to the mea- 
sured time series y az azs and y e /e/j* 

The approximation errors and the gain factors are pre- 
sented in Tables 2-5. Notice that the errors and factors 
are relatively close within each run of the experiment. The 
percentages for each run are of the same order and are 
small. 


V. Conclusions 

This article presents the use of experimental data col- 
lected at the DSS-13 beam waveguide antenna to adjust 
an existing analytical model and to identify a set of single- 
input, single-output models. Four models (azimuth-to- 
azimuth, azimuth- to-elevation, elevation-to-elevation, and 
elevation-to-azimuth) were obtained from the experiment. 
The measurement and analytical techniques used to arrive 
at both of these models are described. The experiment de- 
signed to gather the field data is also described in detail. 

The next step is to obtain a viable and reasonably com- 
pact model of the DSS-13 antenna control system. The 
reduced-order models of the azimuth and elevation axes 
can and will be developed separately. Both models will 
then be used in the design of a model-based antenna po- 
sition loop controller. 


46 


Acknowledgments 


The authors would like to thank Ben Parvin for his support; Jeff Mellstrom for 
his valuable contributions; and Leon Alvarez, Jer-Nan Juang (of NASA Langley 
Research Center), and Michael Thorburn for their helpful discussions. 


References 

[1] W. Gawronski and J. A. Mellstrom, “Modeling and Simulations of the DSS- 
13 Antenna Control System,” The Telecommunications and Data Acquisition 
Progress Report 4 2-1 06, vol April- June 1991 , Jet Propulsion Laboratory, Pasa- 
dena, California, pp. 205-248, August 15, 1991. 

[2] W. Gawronski and J. A. Mellstrom, “Control and Dynamics of the Deep Space 
Network Antennas,” Control and Dynamic System , vol. 57, edited by C. T. Leon- 
des, San Diego, California: Academic Press, 1993. 

[3] J. S. Bendat and A. G. Piersol, Random Data, Analysis and Measurement Pro - 
cedure , 2nd ed., New York: Wiley, 1986. 

[4] J.-N. Juang, L. G. Horta, and M. Phan, “System/Observer/Cont roller Iden- 
tification Toolbox (SOCIT) ,” NASA Technical Memorandum 107566, Langley 
Research Center, Hampton, Virginia, February 1992. 

[5] J.-N. Juang, M. Phan, L. G. Horta, and R. W. Longman, “Identification of Ob- 
server and Kalman Filter Markov Parameters: Theory and Experiment,” Proc, of 
the AIAA Guidance, Navigation , and Control Conf. , New Orleans, Louisiana, 
pp. 1172-1179, August 1991. 


Table 1. Natural frequencies, t Hz, of the four 
transfer functions. 


faxax 

Jaxtl 

Jelcl 

felax 

1.65 

1.62 

2.01 

1.98 

3.24 

3.26 

3.15 

3.21 

4.24 

4.21 

5.22 

5.30 

- 

5.76 

- 

- 


Table 2. 

Approximation errors for the azimuth-to-azimuth 
time series. 

Data set 

Cozoi. percent 

6 ai , percent 

1 

0.31 

0.69 

2 

0.21 

1.00 

3 

0.27 

1.01 

4 

0.27 

1.32 


Table 3. 

Approximation errors for the elevatlon-to-elevatlon 
time series. 

Data set 

telely percent 

6*1, percent 

1 

0.17 

1.07 

2 

0.17 

1.22 

3 

0.23 

2.18 

4 

0.17 

2.05 


Table 4. Approximation errors for the azimuth-to-elevation 
time series. 


Data set 

Cozelt percent 

V 

1 

6.81 

0.0014 

2 

6.75 

0.0013 

3 

6.36 

0.0014 

4 

5,95 

0.0014 


Table 5. Approximation errors for the elevatlon-to-azlmuth 
time series. 

Dataset c«l*z, percent 

w ozb» V 

1 4.12 

0.0014 

2 4.57 

0.0014 

3 4.66 

0.0015 

4 4.19 

0.0014 
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Fig. 1. The experimental setup at the DSS-13 beam waveguide antenna. 







AZIMUTH ENCODER/ELEVATION RATE ELEVATION ENCOOER/AZIMUTH RATE AZIMUTH ENCODER/AZIMUTH RATE 





Fig. 2. Magnitude and phase averages: (a) azimuth rate Input to azimuth encoder output, (b) azimuth rate input to elevation encoder 
output, (c) elevation rate input to azimuth encoder output, (d) elevation rate Input to elevation encoder output, (e) azimuth rate Input to 
azimuth position output phase, and (f) elevation rate Input to elevation position output phase. 
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FREQUENCY, Hz FREQUENCY, Hz 


Fig. 3. Azimuth rata Input to azimuth position output: (a) esti- Fig. 5. Elevation rate Input to azimuth position output: (a) esti- 
mated transfer function from SISO data, (b) system identification mated transfer function from SISO data, (b) system Identification 

model transfer function, and (c) coherence of SISO data. model transfer function, and (c) coherence of SISO data. 



FREQUENCY, Hz FREQUENCY, Hz 


Fig. 4. Azimuth rate Input to elevation position output: (a) esti- Fig. 6. Elevation rate Input to elevation position output: (a) esti- 
mated transfer function from SISO data, (b) system Identification mated transfer function from SISO data, (b) system identification 

model transfer function, and (c) coherence of SISO data. model transfer function, and (c) coherence of SISO data. 
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ELEVATION ENCODER/AZIMUTH RATE AZIMUTH ENCODER/AZIMUTH RATE 








ELEVATION POSITION, deg x 10-3 AZIMUTH 



0 5 10 15 20 25 30 35 


TIME, sec 

|. 8. Detrended azimuth-to-elevatlon time series output data: 
(a) measured and (b) simulated. 



0 5 10 15 20 25 30 35 


TIME, sec 

Fig. 9. Detrended azlmuth-to-azimuth cross-coupling time series 
output data: (a) measured and (b) simulated. 
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Fig. 10. Detrended elevation-to-azlmuth cross-coupling time 
series output data: (a) measured and (b) simulated. 
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Fig. 11. Detrended elevation-to-elevation time series output 
data: (a) measured and (b) simulated. 
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